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ABSTRACT 

We have created a map of the large-scale infrared surface brightness in excess 
of that associated with the atomic interstellar medium, using region-by-region 
correlations between the far-infrared and 21-cm line surface brightness. Our 
study updates and extends a previous attempt with the Infrared Astronomical 
Satellite and Berkeley/Parkes H I surveys. The far-infrared observations used 
here are from the Cosmic Background Explorer Diffuse Infrared Background 
Experiment, which extends far-infrared wavelength coverage to 240 jum, so that 
we are reliably sampling the emission of large, thermal-equilibrium grains that 
dominate the dust mass. The H I data are from the combined Leiden-Dwingeloo 
and Parkes 21-cm line surveys. Using the maps of excess infrared emission at 
100, 140, and 240 /im, we created an atlas and identified the coherent structures. 
These infrared excess clouds can be caused both by dust that is warmer than 
average, or by dust associated with gas other than the atomic interstellar 
medium. We find very few warm clouds, such as the H II region around the 
high-latitude B-type star a Vir. The majority of the infrared excess clouds are 
colder than the average atomic interstellar medium. These clouds are peaks of 
column density, and their excess infrared emission is due to dust associated with 
molecular gas. We identify essentially all known high-latitude molecular clouds 
in the infrared excess maps, and further identify a sample of new clouds with 
similar infrared properties. The infrared excess was correlated with CO line 
brightness, allowing us to measure the ratio of H 2 column density to CO line 
integral (i.e. the N(H. 2 )/W(CO) conversion factor) for high-latitude clouds. The 
atlas of infrared excess clouds may be useful as a guide to regions of relatively 
higher interstellar column density, which might cause high extinction toward 
extragalactic objects at optical and ultraviolet wavelengths and confusion toward 



- 3- 



structures in the cosmic background at infrared and microwave frequencies. 

Subject headings: cosmology: diffuse radiation — infrared: ISM: continuum - 
ISM: clouds — ISM: molecules 
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1. Introduction 

It has been known for some time that interstellar molecular gas exists away from the 
galactic midplane, based on the presence of absorption lines of H 2 and other molecules in 
the spectra of high-latitude stars flSavage et al. 19771 ) an d the presence of millimeter-wave 



emission of CO from regions of high optical extinction (Blitz, Magnani, & Mundy 1984). 
However, the pervasiveness of the molecular component of the local interstellar medium has 
not been convincingly assessed, because of the difficulties of observing molecular gas over 
large areas. Absorption line observations are restricted either to trace elements or to lines of 
sight toward the few hot, bright stars at high latitude, because the electronic ground-state 
transitions of both atomic and molecular hydrogen are in the ultraviolet. Millimeter- wave 
observations of CO are restricted to relatively small areas because all-sky observations are 
at present prohibitively expensive. An early survey of randomly- selected lines of sight found 
a surface filling fraction of 0.5% for CO(l — > 0) emission (Magnani, Lada, & Blitz 1986). 
A large — but, necessarily, incompletely-sampled — survey was recently performed with a 
very low detection rate of 0.3%, suggesting that the northern galactic hemisphere is largely 
'devoid' of molecular gas (Hartmann, Magnani, & Thaddeus 1997). On the other hand, 
there remains the question of whether a survey in a particular observable quantity, such 
as the brightness of collisionally-excited rotational line emission of a particular molecule, 
is really tracing all of the molecular gas. The H2 and CO might have somewhat different 
spatial distributions due to their different photodissociation cross-sections and chemistry, 
and the rotational energy levels of CO may not be excited in low-density (n(H 2 ) < 10 3 
cm -3 ) environments. 

With the advent of the infrared all-sky survey by the Infrared Astronomical Satellite 
(IRAS), a new light was cast on study of the interstellar medium ( [Neugcbaucr et aT 



|1984| , [Low et al. 1984|) . The 100 fim surface brightness was found to be well-correlated 
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with the H I column density on both small and large scales ( Boulangcr and Perault 1988 ) 



demonstrating its value as a tracer of interstellar gas. If the infrared emission arises from 
both the atomic and molecular phases of the interstellar medium, then regions with an 
excess infrared emission relative to the H I column density are likely locations of molecular 
gas. This effect has been demonstrated in Ursa Major (|de Vries et al. 1987|) and for isolated 



cirrus clouds (Heiles, Reach, & Koo 1988; Reach, Koo, & Heiles 1994); in both cases the 
infrared excess was found to be associated with CO emission. Because the infrared surface 
brightness and H I column density are known over the entire sky, it is possible to use the 
difference between the infrared map and an appropriately-scaled map of H I column density 
to produce an all-sky survey of molecular gas. This idea has been exploited by Desert et al. 
(1988), who used the IRAS 100 /zm data and the Berkeley H I survey QHeiles and Habing 
|1974|) to create a catalog of infrared excess clouds. 



In the present paper, we study the distribution and nature of infrared excess clouds 
using relatively recent data from the Cosmic Background Explorer^ (COBE) mission 
( JBoggcss et al. 1992| ) and the Leiden-Dwingeloo H I survey flHartmann and Burton 1997|) . 
These new surveys provide higher sensitivity and higher reliability than the previous 
infrared and H I observations for large-scale emission. More importantly, the COBE 
observations at 100, 140 and 240 /im wavelength are a reliable measure of the emission 
from large, thermal-equilibrium grains that dominate the dust mass. From a detailed study 
of the infrared emission in the Orion region, it was shown that the 100, 140, and 240 fim 



1 The National Aeronautics and Space Administration/ Goddard Space Flight Center 
(NASA/GSFC) is responsible for the design, development, and operation of the Cosmic 
Background Explorer (COBE). Scientific guidance is provided by the COBE Science Working 
Group. GSFC is also responsible for the development of the analysis software and for the 
production of the mission data sets. 
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emission sample essentially the same dust temperature along the line of sight ( Wall et"aT 



|1996|) . Our first results ( |Reach et al. 1994]) , based on comparing the COBE 240 /im optical 



depth to the Berkeley H I surveys, encouraged us to pursue a more thorough comparison of 
the H I and infrared data. 

The large-scale distribution of molecular gas is important for a number of practical 
applications. 

Extinction of extragalactic objects — Extragalactic observations at visible and shorter 
wavelengths are affected by extinction even at high galactic latitude. In order to estimate 
the extinction, which affects both the brightness and color of extragalactic objects, it 
has been necessary to rely on 21-cm line surveys (cf. Purstein and Heiles 1978|) . If an 
extragalactic object lies behind a high-latitude molecular cloud, its extinction will be 
significantly underestimated using the 21-cm line surveys. The results derived here will 
help observers to identify regions where anomalously high extinction from molecular clouds 
can be expected. A recent effort by another group (Schlegel, Finkbeiner, & Davis 1997) 
systematically addresses this issue, using the COBE and IRAS data to create a map of the 
extinction. Large-scale variations in the temperature and gas-to-dust ratio were calibrated 
by Schlegel et al. in a manner very similar to ours. By calculating the column density using 
large-scale average dust temperature, it is likely that the Schlegel et al. dust maps will 
underestimate the extinction toward relatively cold clouds such as the molecular clouds we 
have identified. We therefore recommend that the Schlegel et al. maps be used to estimate 
the extinction, and that our molecular cloud atlas be used as a supplement, warning of 
cold, high-extinction clouds. Lines of sight behind these clouds should certainly be avoided 
unless 1-3 magnitudes (Ay) of extinction can be tolerated. 

Shadowing of distant X-rays — Interstellar clouds produce distinct shadows on the soft 
X-ray emission from the Galactic halo and the extragalactic emission ( [McCammon and 
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Sanders 1990). These shadows have been compared to 21-cm line maps and to IRAS 100 
fim surface brightness maps (Snowden, McCammon, & Verter 1993; [Wang and Yu 1995| ). 



In both cases, the interstellar column density is underestimated in the presence of molecular 
clouds, because molecular gas can absorb X-rays and molecular clouds are relatively faint 
at 100 /j,m (and nearly invisible at 60 /mi; [Laureijs et al. 1996j ), compared to atomic gas. 



We provide both a map of infrared excess clouds and a calibration of their column density 
in the present paper. 

Diffuse 'y-ray emission and the N(H 2 )/W(CO) factor — A significant source of the 7-ray 
surface brightness of the sky is due to interaction of cosmic rays with interstellar gas. The 
expected near-linear proportionality of 7-rays with total column density allows a calibration 
of the molecular column density, when correlated with H I and CO maps ([Strong et al. 



1988| ; |Digel et al. 1996| ). In the presence of H 2 clouds with relatively little CO emission, 



this calibration can be significantly biased. Here we calibrate the molecular column density 
using the infrared excess, and compare the results to 7-ray studies. The ratio of H2 column 
density to CO line integral is important for assessing the vertical mass distribution of 
molecular gas in the disk of our Galaxy. 

Relation to external galaxies — Low-metallicity galaxies such as the Large and Small 
Magellanic clouds contain relatively fewer giant molecular clouds, and more translucent 
regions — due to the lack of dust. High-latitude molecular clouds such as studied here share 
empirical similarities with low-metallicity galaxies, and should provide insight into the 
nature of the latter. 



2. The Correlation of Far-Infrared Emission with H I 
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2.1. Observations 

The far-infrared data used in this work originate from the Diffuse Infrared Background 
Experiment (DIRBE, [Hauser et al. 1997]) . The DIRBE produced a redundantly sampled 



all-sky survey in 10 wavebands from the near to far-infrared; here we concentrate on the 100 
/jiii, 140 (Urn, and 240 /im wavebands, which are dominated by the emission from interstellar 
dust. In particular, these wavebands are dominated by the emission from the larger grains, 
in thermal equilibrium with the interstellar radiation field, that dominate the mass of dust 
in the interstellar medium. Nonetheless, a significant source of large-scale emission, even 
in the far- infrared, is from the zodiacal light. For this work, we use data from which a 
model for the zodiacal light has been subtracted ( [Reach et al. 1996| ; [Kelsall et al. 199~8| ). 



To minimize detector noise, which is significant for the 140 /im and 240 /mi wavebands, 
we use the average of the zodiacal-light-subtracted weekly maps over the cold mission (44 
weeks). The instantaneous beamsize of DIRBE is 42'. In order to avoid confusion between 
the interstellar clouds that we hope to study and bright infrared point sources, we use 
data from which all pixels containing bright and apparently unresolved structure, as well 
as circular regions around the Magellanic clouds, have been masked ( |Arendt et al. 191H5 ). 



Furthermore, we will not consider here the portion of the sky close to the galactic plane 
(\b\ < 20°), where clouds from a wide range of distances overlap on the sky. The DIRBE 
data are gridded into a sky-cube projection, which we degraded to a pixel size of 39' so 
that each pixel is nearly independent. We assume that the point-to-point uncertainty 
in the infrared surface brightness at (100, 140, 240 /an) is a quadrature combination of 
(0.15, 1.3, 0.75 MJy sr~ x ) and 2% of the surface brightness, to account for detector noise, 
pixelization uncertainties, source confusion, and gain drifts. These uncertainty estimates 
were determined from inspection of the data. At 140 and 240 /im the point-to-point 
uncertainty is comparable to the expected instrument noise (cf. [Hauser et al. 199~8| ), but at 
100 /xm it is significantly larger, so we must have included some of the detected small-scale 
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sky structure of the interstellar medium. The uncertainties are important in the present 
work mostly for weighting the data, but they do scale the \ 2 statistics. 

The H I 21-cm line data used in this work originate from the Leiden-Dwingeloo survey 
of the Northern sky ( |Hartmann and Burton 1997|) and the Parkes survey of the Southern 
sky ( |Cleary et al. 1979|) . The Dwingeloo 25-m telescope has a main beam size of 36', which 
is well-matched to the DIRBE beamsize. An important aspect of the Leiden-Dwingeloo 
survey is the removal of stray radiation, which is up to half of the observed signal at high 
galactic latitudes. The near and far sidelobes of the telescope response were mapped using 
bright point sources, then the sidelobe emission was estimated by convolving the sky map 
with the beam ( Hartmann et al. 1996| ). For this work we use the integrated 21-cm line 



brightness, converted into H I column density assuming the emission is optically thin. This 
latter assumption is unlikely to be violated at high galactic latitude unless the H I gas 
is colder than 20 K. Gas this cold must be well-shielded from the interstellar radiation 
field, or else photoelectric heating will raise its temperature to at least 50 K at typical 
interstellar pressure. In such cold, shielded regions, the H is rapidly converted into H2, so 
that the H I column density becomes a negligible fraction of the total. We integrated the 
H I line over the range of -100 to +100 km s _1 , which includes the bulk of local interstellar 
gas but excludes high-velocity clouds. The high-velocity clouds have been shown to be 
deficient in interstellar dust ( |Wakker and Boulanger 1986Q ; therefore they can be ignored 
in the present study. For declinations south of —30°, we used the 21-cm survey performed 
with the Parkes 64-m telescope flCleary et al. 1979|) . The Parkes survey provides critical 



information on the southern sky, but it does suffer from undersampling (every beamwidth in 
right ascension but every other beamwidth in declination) and incompleteness (some survey 
scans missing). Furthermore, the Parkes survey has not been corrected for stray radiation. 
We checked the calibration of the Parkes and Leiden-Dwingeloo surveys by comparing their 
column densities where the surveys overlap; the calibrations were in good agreement. We 
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assume the point-to-point uncertainty in the H I column density is a quadrature sum of 
1 x 10 19 cm -2 and 2% of the column density, to account both for uncertainties in the stray 
radiation and the gain. The uncertainties were not rigorously determined; they are based 
on the estimated gain and stray radiation accuracy (|Hartmann and Burton 1997|) . 

2.2. Region- by-region correlations 

EDITOR: PLACE FIGURE HERE. 

In order to identify molecular clouds in the infrared maps, we must first remove the 
infrared emission from dust associated with the atomic interstellar medium. We presume the 
column density of the atomic interstellar medium is linearly traced by the integrated 21-cm 
line brightness. The infrared emission, however, is not a linear tracer of the interstellar dust 
because it is sensitive to the dust temperature. We considered two methods to calibrate 
the dust temperature variations: (1) calculate the column density at each pixel using the 
100, 140, and 240 /im brightness ratios assuming a single dust temperature along the line 
of sight, and (2) calculate the ratio of infrared to H I surface brightness from region to 
region on the sky. The first method is simpler and, in principle, it could detect temperature 
changes on angular scales as fine as 1°. But in practice the 140 and 240 /im maps are 
very noisy and the 100 /zm map suffers from large-angular-scale zodiacal light residuals 
(and all wavelengths contain some cosmic infrared background radiation). Furthermore, we 
show below that the dust associated with molecular gas is significantly colder than that 
associated with atomic gas, so that a single dust temperature does not apply to the lines 
of sight we are most interested in. Therefore we used method (2), by dividing the sky into 
'cells', with 10° radius, on a regular grid every 10°. Within each cell, we measured the slope 
and intercept of the pixel-to-pixel infrared surface brightness as a function of H I column 
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density. To illustrate the procedure, Figure [1] shows an example. In the upper panel, 
we show the 240 /im surface brightness versus H I column density for the entire dynamic 
range of both quantities, within the region 15° in radius centered on galactic coordinates 
(/,&) = (140°, 35°). At low H I column densities, there is a linear rise of infrared surface 
brightness with H I column densities, while at higher column densities there is a wide 
dispersion. This dispersion is due to infrared emission associated with gas other than H I. 
At low column densities, interstellar H remains atomic, but where the column density is 
large, the H becomes molecular. For a uniform- density cloud, the H2 column density is 
predicted to exceed the H I column density when 

N(E I) > 2.5 x 10 20 ( -) 2 ( 1 cm" 2 (1) 

v ; V100 cm" 3 / V80 KJ v ; 

( [Reach et al. 1994] ). In low-density regions, H2 cannot form, while in denser clouds such as 

account for the 21-cm absorption features in front of radio sources (see review by |Kulkarni 



and Heiles 1988|) , the gas readily forms molecules at column densities > 3 x 10 20 cm -2 . 
Figure [l] shows this effect: the points at higher H I column densities deviate systematically 
upward from the straight line fitted to the lower column-density points. This upward 
deviation, presumably due to infrared emission from dust associated with molecular gas, is 
what we will call the 'infrared excess.' We presume that on large scales, the warm and cold 
phases of the interstellar atomic gas exist in pressure equilibrium, and the cold clouds have 
a relatively small filling factor, so that all lines of sight cross warm regions and some (or 
all) cross cold regions. In order to exclude cold clouds with column densities sufficient to 
form H2 from our infrared-H I correlations — which are intended to measure the properties 
of the atomic medium only — we exclude all points with a column density more than 
3 x 10 20 cm -2 above the minimum column density for the 'cell.' The minimum column 
density is presumed to represent the warm H I, and by allowing a different threshold in each 
cell we allow for the longer path length through the warm gas Galaxy at low latitudes. The 
threshold column density for the Ursa Major/Ursa Minor/ Camelopardalis region is shown 
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vertical dashed line in Figure [l|. 

EDITOR: PLACE TABLE g HERE. 

For the low-column density pixels within a given sky 'cell,' we perform a linear fit of 
infrared surface brightness versus H I column density, taking into account the uncertainties 
in both quantities. For the sample cell, the points used in the fit and their uncertainties are 
shown in the lower panel of Figure |I|. Even at low column densities, there are some outlying 
positive points, but in no case do the fits appear to be biased by outliers. The resulting 
slope and intercept for two sample cells, widely separated on the sky, are shown in Table |I[ 
In both cases the fits are 'good' in the sense that the reduced xt is essentially unity (to 
within our understanding of the measurement uncertainties). At 240 /im, the infrared-H I 
slopes of the two cells are similar, but at 100 /zm the slopes are significantly different. Using 
the ratio of slopes at 100 and 240 /im, (and assuming the interstellar dust emissivity varies 
as v 2 ; Draine and Lee 1984 ), we find that the temperature of the dust associated with H I 



is 17.2 ± 0.2 K in the cell centered on (140°, 35°), and 19.0 ± 0.3 K in the cell centered 
on (40°, 35°). Even though a 2 K temperature difference between the two cells may seem 
unimportant, a change in temperature from 17 to 19 K increases the emission at 100 /im 
by a factor of 2.4 and the 240 /zm emissivity by a factor of 1.7, so that the ratio of 100/240 
/j,m emission increases by 60%. The ratio of dust optical depth at 100 fim (assuming u 2 
emissivity) to gas column density in the two cells are 6.4 and 4.6, respectively, in units 
of 10~ 25 cm 2 . The absolute calibration uncertainty is ±15% in these dust-to-gas ratios, 
but the uncertainty in cell-to-cell variations is much less. Therefore, the (30%) difference 
between the dust-to-gas ratios in these cells is statistically significant. 



EDITOR: PLACE FIGURE | HERE. 
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We have calibrated the differences in the dust temperature and dust-to-gas ratio 
from region to region by performing independent fits of the infrared surface brightness 
as a function of H I column density at each wavelength and in each sky cell. These fit 
parameters were well-constrained for most of the sky, but there were some regions where 
the infrared emission was relatively poorly-correlated with the H I, even at low column 
densities. If a molecular cloud complex or dust hot-spot is comparable to or larger than 
our cell size (10°), then it is not possible to calibrate the infrared-H I slope. These regions 
are generally near the galactic plane, and for the most part they can be identified with 
well-known star-forming complexes. For each of these complexes, a detailed study using 
higher-resolution maps is needed to sort out the variations in gas phase and dust heating 
due to embedded stars (e.g. |Boulanger et al. 1998| for the Chameleon complex). 



All-sky maps of the infrared-H I slope and offset at 100 /mi are shown in Figure 0. In 
cells where the fits were found to be poor, we interpolated using adjacent cells. In cells 
where the fit is good, we interpolated to a finer, 37' grid in order to avoid discontinuities at 
cell boundaries. Two regions that stand out in the infrared-H I slope are at (300°, +50°) 
and (90°, —40°). The former is due to the effect of the nearby, early-type star Spica, whose 
situation at high latitudes allows it to heat dust over a very large apparent area ( |RcynolcTs 



55]), while the latter is due to the MBM 53/54/55 complex of molecular clouds ( |Magnani et| 



al. 1985| ). Our cell size was chosen so that regions like these would not completely disappear 



once the infrared emission associated with H I is subtracted. In addition to these relatively 
compact regions, there are real variations in the 100 micron-HI slope due to temperature 
variations from place to place of the interstellar dust associated with atomic gas. These 
variations appear in Fig. |2] as large-scale structures. One noticeable problem is at southern 
declinations, where the H I data from the Parkes survey was used. The infrared emission 
per H atom is lower there than average, which could be due in part to stray radiation effects 
in the Parkes survey. 
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2.3. Implications for the far-infrared background 

EDITOR: PLACE FIGURE 131 HERE. 



EDITOR: PLACE FIGURE | HERE. 

For each of the far-infrared wavelengths we have considered here, we find that a 
positive residual sky brightness remains after subtracting the infrared emission correlated 
with the H I column densities. Histograms of all the zero-intercepts of the infrared-H I 
correlation for out grid of sky cells are shown in Fig. |j; it is clear that the zero-intercept is 
consistently positive. This same result has also been found from the correlations of infrared 
brightness with H I column density in the Lockman Hole and around the North Ecliptic 
Pole ([Arendt et al. 19981) , and for the entire low-column-density portion of the northern sky 



(|Boulanger et al. 1996[) . At least part of the residual emission is related to the Solar System, 



because one can see the distinctive pattern of the ecliptic plane in Fig. [| However, the 
amplitude of the residual zodiacal light, based on its ecliptic latitude dependence, is small. 
Furthermore, the spectrum is inconsistent with the zodiacal light, being relatively bright at 
240 /mi, which argues against a Solar System origin ( Pwek et al. 1998|) . The far- infrared 
emission remaining after the subtraction of the emission from interstellar dust could be the 
extragalactic background due to unresolved galaxies, with obvious cosmological significance; 
alternatively, it could be due to a component of our galaxy with infrared emission but little 
atomic gas. We will address the cosmological issue only briefly here, as it is the primary 
focus of other recent papers ( [Puget et al. 1995| ; |Hauser et al. 199"S| ). In the present work, 



we have allowed for the possibility that the infrared emission per unit H I column density 
varies from place to place. In the presence of uncertainties in both variables, the slope and 
zero-intercept of a correlation tend to be correlated. Because the H I column density is 
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always positive, the slope and intercept are negatively correlated; furthermore, the intercept 
is biased toward positive values where the correlation is weak. Comparing the slope and 
offset of all of the sky cells with a significant fit, we find a near-perfect anti-correlation. 
This could potentially lead to a false detection of the cosmological background, especially 
if, at low column densities, the interstellar dust is associated with an interstellar medium 
other than atomic gas. 

EDITOR: PLACE TABLE @ HERE. 

The infrared emission per unit H I column density has a well-defined mean value, with 
few regions of abnormally high or low values. Table [| lists the mean and rms dispersion of 
the infrared emission per unit H I column density and the zero-intercept of the infrared-H I 
correlation at 100, 140, and 240 fim. It is clear that the range of slopes and intercepts 
obtained over the high-latitude sky is significantly larger than the uncertainty of an 
individual value, so that the range of slopes is not due to random measurement or fitting 
uncertainties. The wavelength-dependence of the slopes in the two sample regions (Tab. |TJ) 
indicates that this is mostly a dust temperature effect. This is a gradual, large-scale effect, 
that is not confined to the few outstanding cells. The cells with an infrared-H I slope more 
than 1.5 times the rms away from the mean are shaded in Fig. [3] and Fig. |j. It is evident 
that the outlying regions in terms of infrared-H I slope are also outlying in terms of the 
zero-intercept (as expected because of the anti-correlation of slope and intercept mentioned 
above). Over the entire sky, there is a significant, positive sky brightness that cannot be 
explained by interstellar dust associated with H I. If the range of slopes shown in Fig. [3] is a 
real effect due to the interstellar medium, then the range of offsets in Fig. |] is also due to 
the Galaxy. The inferred brightness of an isotropic extragalactic background must then be 
lower than the lowest offset value. In principle, variations in the infrared-H I correlation on 
scales smaller than those chosen for the present work could lead to an even wider range of 
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slopes and offsets. The dispersion of the putative infrared background brightnesses, listed 
in the last column of Table [| precludes more than a '2-er' (loosely speaking) confidence level 
detection of the infrared background, unless more restrictive assumptions about the galactic 
emission are made. The significance of this effect is discussed in more detail by Hauser et al. 
(1998) and Shafer et al. (1998). Fortunately, the separation of the galactic and cosmological 
contributions to the far-infrared sky brightness can be resolved using higher-resolution and 
sensitivity far-infrared observations, such as have recently been performed with the Infrared 
Space Observatory ( Kawara et al. 1997| , [Puget et al. 1998] ) and will eventually become 
possible with the NASA Space Infrared Telescope Facility and the ESA Planck Surveyor. 
These observations should reveal whether the brightness of the faintest parts of the sky is 
produced by an ensemble of extragalactic sources or by interstellar cirrus. 



3. The Distribution of Infrared Excess 

3.1. Definition of infrared excess 

The infrared excess is defined to be the observed infrared surface brightness minus the 
contributions from dust associated with the atomic interstellar medium, the zodiacal light, 
and the cosmic infrared background. We calculate the infrared excess using 

I? = h - S * N(R I) - O, (2) 

where S and O are the locally-interpolated slope and offset of the infrared-H I correlation. 
The offset O contains residual zodiacal light and the cosmic infrared background, and is 
not useful for measuring properties of the interstellar medium. We created infrared excess 
maps independently for the 100, 140, and 240 fim wavebands, so that we can measure the 
far-infrared colors of the infrared excess. The infrared excess map at 240 fim is shown in 
an all-sky projection in Figure |5]. The infrared excess is 'patchy' with significant structure 
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on scales smaller than the size of the 'cells' within which the infrared-H I slopes were 
determined. We have checked that the maps are not very different if the cell size is twice 
as large as the one we finally adopted. The significance of a given value of the infrared 
excess can be determined statistically from its probability distribution. Figure | superposes 
the histograms of the positive and negative pixels in the infrared excess map. It is clear 
that the probability distribution is not random, and that the positive values significantly 
outnumber the negative values when /fgo > 0-3 MJy sr -1 , which we will use as the lower 
limit to define infrared excess clouds. 

EDITOR: PLACE FIGURE | HERE. 
EDITOR: PLACE FIGURE | HERE. 

3.2. Atlas of infrared excess clouds 

A table of clouds would not suffice to describe the infrared excess, because the emission 
has some continuity over large scales. High-latitude molecular clouds have been shown 
to be associated with large H I features (Gir, Blitz, & Magnani 1994), suggesting their 
formation is related to events or processes (like supernovae and stellar winds) that organize 
the interstellar medium on large scales. For presentation here, and practical utility to 
other researchers, we present the infrared excess as an atlas consisting of 8 maps with 
simple coordinate projections. We have preferred a Cartesian projection for latitudes 
20° < |6| < 50°, so that positions can be easily located with a ruler (Fig. [7]). We also made 
maps of the two galactic polar caps, in a simple orthographic projection (Fig. The 
maps we present here are for 100 /im wavelength, the highest- sensitivity of the DIRBE 
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far-infrared wavebands. The maps at 140 and 240 /im are generally similar. At 100 /im, 
there are more infrared excess clouds, because both warmer-than- average and molecular 
clouds produce infrared excess at 100 /im. However, the vast majority of structures in the 
far-infrared excess maps are shown in the next sections to be cold. 

EDITOR: PLACE FIGURE HERE. 

EDITOR: PLACE FIGURE § HERE. 

In order to study the properties of the 'clouds' and regions of sky with infrared excess, 
we have identified two lists of positions and apertures within which to perform photometry. 
The first list is the catalog of previously-known molecular clouds compiled by Magnani, 
Hartmann, & Speck (1997). (We combined entries that would be redundant at our low 
resolution.) These are indicated on the individual panels of the atlas. Second, we made a 
list of some of the remaining peaks in the infrared excess maps, with brightness comparable 
to the known molecular clouds. This list was created by hand, upon inspection of both 100 
/im and 240 /im infrared excess maps, and it is incomplete. The 'new' clouds were named 
in the format DlKUl ± bb where III is the galactic longitude and bb is the galactic latitude; 
the prefix 'DIR' stands for 'diffuse infrared' as in the name DIRBE. Third, for reference, we 
labeled many peaks in the infrared excess maps due to bright point sources, most of which 
are nearby galaxies (Odenwald, Newmark, & Smoot 1997). The list of previously-known 
clouds is shown in Table |^, and the newly- identified clouds are shown in Table |j. For each 
cloud, we averaged the infrared excess surface brightness at 100 and 240 /im within the 
cloud boundary. Only clouds detectable both at 100 /im and 240 /im are listed. The total 
infrared excess flux at 100 /im, the color temperature obtained from the 100/240 /im ratio 
(assuming emissivity proportional to z/ 2 ), and the angular diameter of each cloud is shown 
in Tables § and £|. 
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EDITOR: PLACE TABLE | HERE. 



EDITOR: PLACE TABLE | HERE. 

The infrared excess clouds in general are relatively bright at 240 /iin, and their 
far-infrared color temperature is significantly lower than that of average atomic clouds. This 



is illustrated in Figure [10], where histograms of the color temperatures of infrared clouds 
are shown. The fact that the vast majority of infrared excess clouds are cold confirms that 
the infrared excess is generally not due to clouds that are warmer than average. Indeed, 
because the clouds are colder than average, they tend to be under-luminous in infrared 
emission. This effect has also been found by Lagache et al. (1997), who showed that 
molecular regions at high latitude can be identified by (1) deficiency in 60 /zm emission 
and (2) low 140/240 /im color temperatures relative to atomic cirrus. One of the main 
uncertainties that has been faced by past studies of infrared excess as a tracer of molecular 
gas (e.g. |Desert et al. 1988 ; [Reach et al. 199~4f) was the uncertainty as to whether the 



infrared surface brightness at 60 /xm, observed by IRAS and often compared to the 100 
jum brightness, traces variations in column density or dust temperature. Using the COBE 
observations at 240 /xm, this problem has now been resolved. 



3.3. Warm infrared excess clouds 

While they are rare, there are some infrared excess clouds that are warmer than 
average. The single most prominent example is the large H II region around the nearby 
B star a Vir, or Spica, which was first discovered by its extensive Ha emission ( |Reynolds| 
|S5J). The Spica H II region is exceptional, in that it is produced by one of the closest — at 
d = 80 ± 5 pc (ESA 1997|) — early-type stars; it is at high galactic latitude, where it can be 
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easily separated from other sources; and the star is presently in a region of relatively low 
interstellar medium density, so that its H II region covers a region some 15° across. The 
dust in the Spica H II region is significantly warmer than typical interstellar cirrus, and 
it is unusual in that it also has 12 and 25 //m emission significantly higher than expected 
from the 100-240 /zm emission, suggesting a modified grain size distribution flZagury 19971 ). 
There are no other 'warm' infrared excess clouds comparable to the Spica H II region in size 
and brightness, most likely because no other early-type star is as advantageously situated. 

One significant 'warm' infrared excess cloud that has not been noted before is 
DIR015+54. This is one of the brightest 100 /im infrared excess clouds at high positive 
galactic latitudes (Fig. |8|a), but it is not detected at 240 /im (and therefore is absent 
from Table |j). In order to see the structure of this cloud at higher angular resolution, 
we extracted IRAS images at 12, 25, 60, and 100 /im from the IRAS Sky Survey Atlas 
(ISSA). The cloud is relatively bright at 60 /im, with a local-background-subtracted 
surface brightness, averaged over the central 30', of 0.8 MJy sr" 1 . The ratio of 60/100 
/im brightnesses is 0.46, a factor of 2 higher than that of atomic 'cirrus' and much higher 
than that of molecular clouds. The ISSA 60 /im image is shown in Figure |||. The cloud is 
structured in a form qualitatively suggestive of a bow-shock produced by an object moving 
toward the northeast. Bow-shocks were discovered in the IRAS data due to the stand-off 
between the interstellar medium and the stellar winds of early-type stars with significant 
transverse velocities ( |van Buren and McCray 1988|) , and their morphologies and 60/100 
/im colors are comparable to the cloud shown in Figure |9|. However, there is no early-type 
star in the vicinity that could plausibly create this nebula, either as an H II region or a 
bow shock. Such a star would be very bright at visible and ultraviolet wavelengths, owing 
to the low extinction along the line of sight and the proximity needed to heat a region of 
such a large angular size. There is a peak of 60 fim emission approximately at the location 
where a wind-source would be located. The nature of this emission peak, listed in the IRAS 
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Faint Source Catalog as F15101+1206 ( |Moshir et al. 1989 ), is presently unknown, though 
statistically it is likely to be extragalactic (Reach, Heiles, & Koo 1993). Apart from some 
of the peaks in the nebula, which appear in the IRAS Faint Source Catalog, nearby there 
are only some galaxies and stars of type A and later in the SIMBAD database. We find 
that there is a faint optical counterpart of F15101+1206 in the Palomar Observatory Sky 
Survey^, consistent with the idea that F15101+1206 is extragalactic. There is no optical 
counterpart for the diffuse emission of DIR015+54. The infrared colors of DIR015+54 are 
similar to those of the nearby galaxy NGC 6822 ( [Israel et al. 1996 ), suggesting perhaps 
that the cloud could in fact be another low surface brightness member of the Local Group. 
We inspected the individual H I spectra from the Leiden-Dwingeloo 21-cm line survey near 
this position, and we find that the cloud is clearly present, with a 'normal' LSR velocity 
of -17 km s -1 and both narrow and wide line components (dispersions 3 and 15 km s" 1 , 
respectively). We suspect therefore that DIR015+54 is a Galactic cloud that happens to 
be somewhat brighter than average in the far-infrared, and much brighter than average at 
60 fim. The nature of this cloud and the other large, 'warm' clouds in the COBE/DIKBE 
and IRAS data remains mysterious. They are either locally-heated regions or regions with 
unusual grain size distributions, but we do not know why; this should be investigated 
further. 



2 Based on photographic data of the National Geographic Society - Palomar Observatory 
Sky Survey (NGS-POSS) obtained using the Oschin Telescope on Palomar Mountain. The 
NGS-POSS was funded by a grant from the National Geographic Society to the California 
Institute of Technology. The plates were processed into the present compressed digital form 
with their permission. The Digitized Sky Survey was produced at the Space Telescope 
Science Institute under US Government grant NAG W-2166. 
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3.4. Comparison with known molecular clouds 

Our atlas of infrared excess clouds reveals essentially all previously known high-latitude 
molecular clouds in the literature flMagnani et al. 1996|) . First, consider the clouds 



originally found by Magnani, Blitz, & Mundy (1985, hereafter MBM); these clouds were 
first identified as regions of extinction on the Palomar Observatory Sky Survey plates, and 
then they were detected in CO (J = 1 — > 0) emission with radio telescopes. At our low 
angular resolution, the MBM clouds comprise essentially 36 individual regions. Of these, 
well-formed clouds with significant infrared excess (i^o > 1 MJy sr _1 ) were clearly detected 
from 34 regions, for a 94% detection rate. MBM 19 and MBM 9, the two regions for which 
we do not see an infrared 'cloud', both are in regions of significant infrared excess, but there 
is no well-defined cloud at their position. Indeed this is somewhat true of some of the other 
clouds as well: the center of the infrared excess region does not coincide with the center 
of the MBM cloud. We suspect that the MBM positions are relatively dense cores (with 
well-defined extinction on the optical plates), while the overall cloud structure is better 
represented by the infrared excess map. 



EDITOR: PLACE FIGURE 10 HERE. 



Our atlas of infrared excess is significantly different from the catalog of 515 infrared 
excess clouds obtained by Desert et al. (1988) [DBB]. We have attempted to associate 
each entry of Tables |3] and |] with a cloud in the DBB catalog. For 62 % of the MBM 
clouds, we find corresponding entries in the DBB catalog. This is somewhat higher than 
the 47% detection rate determined by DBB, because of the way we grouped MBM clouds 
into complexes and considered a match if the DBB and MBM positions fall within the same 
cloud in Figs. |7] or |8|. Specifically, this means we consider the MBM clouds to be much 
larger than the high-extinction and CO-emitting cores. The higher detection rate of MBM 
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clouds in our atlas is due to the higher quality input data that we used; both infrared 
and HI are of substantially higher quality than the data available to DBB. Many of the 
clouds that are not listed in DBB are very bright in our atlas, including both well-known 
(MBM) and new clouds. On the other hand, a significant fraction of DBB clouds are not 
confirmed by the present analysis. In total, we find counterparts for 45% of the DBB clouds 
and hints for another 7%, leaving 48% with no evident infrared excess in our maps. Of 
the unconfirmed clouds, 58% have a 'significance' criterion less than 5 listed in the DBB 
table, but 10% (mostly H I artifacts) are listed with 'significance' greater than 7. The 
uncertainties in the infrared excess method — especially with the lower-quality infrared and 
H I data used by DBB, but also with the results presented in this paper — are dominated 
by systematic uncertainties rather than random noise. We have therefore chosen to make a 
relatively short table of 81 selected clouds, and to present the sky maps for use by others. 



EDITOR: PLACE FIGURE 11 HERE 



The interstellar medium is not an ensemble of clouds with well-defined boundaries — it 
is coherent over large angular scales. Therefore a comparison of cloud catalogs can be 
misleading. There are few regions where the CO(J = 1 — > 0) line has been mapped in 
emission on scales large enough to make a meaningful comparison with our atlas. The 
largest-scale high-latitude CO map of which we are aware was made with a 1.2-m telescope, 
whose large beam (9') is well-suited to comparison with our results. A large portion of 
the Ursa Minor-Ursa Major region has been covered in a number of surveys, including 
one very large, regularly sampled grid ( Heithausen et al. 1993|) . We have projected our 
infrared excess map onto the grid of the CO map assembled by Heithausen et al. (1993) 
for detailed comparison — see Figure |TT|. The agreement between the contour of integrated 
CO line brightness, W(CO), and the infrared excess contour is remarkably good, and 



the infrared excess map contains no significant clouds in regions without CO. Figure 12 
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shows a pixel- by-pixel scatter diagram for these maps. A linear fit, taking into account 
the statistical uncertainties of 0.5 K km s _1 in W(CO) and 0.3 MJy sr" 1 in yields the 
following: 

W(CO) = (-0.19 ± 0.05) + (1.17 ± 0.05) (j^/MJy sr" 1 ) K km s"\ (3) 

where only statistical uncertainties were accounted for. The correlation of 100 /im emission 
with H I column density in this region (which largely overlaps with that shown in Fig. |]) 
yields 

Jioo = 0.47 ± 0.1 + (0.62 ± 0.03) (iV(H I)/10 20 cm~ 2 ) MJy sr _1 . (4) 

If the infrared emission per H nucleus were the same in atomic and molecular gas, then we 
find that the H 2 column density per unit integrated CO(l — > 0) line brightness is 

iV(H 2 )/W(CO) = (0.7 ± 0.1) x 10 20 cm" 2 K" 1 km" 1 s. (5) 

If they were sufficiently orthogonal, we could determine the separate emissivities of 
the atomic and molecular gas using a multiple linear regression between the infrared, H I, 
and CO surface maps over the parts of the region that have been observed at all three 
wavelengths. Table [|] shows the slopes from fits of the form 

I v = aN(R I) + bW(CO) + c (6) 

for several infrared wavelengths. It is evident that a and b depend on wavelength differently, 
and the interpretation is straightforward: the dust mixed with the atomic gas has a different 
temperature from the dust mixed with the molecular gas. The color temperature inferred 
for the dust mixed with the atomic gas is T\ = 17.0 ± 0.3 K, and for the dust mixed with 
the molecular gas it is T 2 = 14.7 ±0.2 K. The observation that atomic and molecular clouds 
have different infrared spectra implies that their relative column densities cannot be simply 
determined using observations at a single wavelength. This has an important practical 
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implication. Naively, the ratio 

X naive = 0.5b/a (7) 

gives the column density of H 2 per unit CO line integral; values of X naive are listed in the 
last column of Table |5|. It is evident that the naively inferred value of X changes by a factor 
of 2 depending on the wavelength, which is inconsistent with the idea that the infrared 
emission is a linear tracer of the total column density. We can, however, calibrate this effect 
using the results obtained here. We assume that within the molecular phase, the gas is 
completely molecular and the dust-to-gas ratio is the same as in the atomic phase. Then 
the column density of H 2 per unit CO line integral is 

( |Magnani and Onello 1995|) which results in the last column of Table ||. When calibrated 
this way, a self-consistent description of the infrared emission from the atomic and molecular 
phases could be obtained from observations at a single wavelength — if the temperature 
contrast between the atomic and molecular components were the same everywhere. Clearly, 
observations at 240 fim are the best infrared tracer of molecular gas, because the 240 /zm 
emission is less sensitive to temperature changes. 

The H 2 column density per unit CO line integrated intensity has been measured by 
several techniques and for several different regions, with a fairly wide range of results. The 
'canonical' method is to correlate 7-ray emission — due to interaction of cosmic rays with 
nuclei — with the CO and H I emission simultaneously. This correlation has been performed 
in the galactic plane, where the CO emission is dominated by giant molecular clouds, and 
X = (2.3 ± 0.3) has been found f\ ( [Strong et al. 19881) ■ The correlation of 7-ray emission 



and CO line brightness in the Ursa Major/Ursa Minor region yielded X = (0.92 ± 0.1) 



5 Units for X throughout this paper are 10 20 cm 2 [K km s 1 1 
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( pigel et al. 1996 ). The 7-ray method, like the 'naive' infrared method, assumes that the 



7-ray emission properties of the atomic and molecular clouds are the same. In fact, recent 
observations of lower-energy 7-rays indicate that the value of X inferred depends on the 
energy of the 7-rays, with values ranging from (1.1-4) flStrong et al. 1994]) . If the difference 
is due to processes that only operate at low energies, then the 'canonical' value of 2.3 holds, 
for it was determined at high energies; however, until this is understood, there remains the 
possibility that even the 'canonical' X is uncertain by a factor of 2. 

The infrared method of locating and estimating the column density of H 2 is especially 
useful because a wide range of physically different types of cloud can be observed. For 
the galactic center, Sodroski et al. (1995) found that the value of X is an order of 
magnitude lower than in the outer Galaxy. Their analysis takes into account the higher 
dust temperature in the galactic center, because they used the dust optical depth from 
140 and 240 //m observations. At high latitudes, the molecular clouds are smaller, so one 
might anticipate a difference in the H2 column density per unit CO line integrated intensity 
due to the higher radiation fields that can penetrate a smaller cloud. An analysis very 
similar to our multiple linear regression was performed for the Ursa Major clouds using 
higher-angular-resolution 21-cm line observations with the Effelsberg 100-m telescope and 
IRAS 100 /iin data ( |de Vries et al. 1987| ). Their infrared map clearly revealed the presence 



of the clouds that were not evident in the 21-cm line map, and the combined H I and CO 
maps were found to reproduce the infrared brightness, with X naive = (0.5 ± 0.3). Using the 
temperature of atomic and molecular clouds we found for the entire North Celestial Pole 
region, 

X = 3.8X naive = (1.9 ± 1.1) x 10 20 cm- 2 [K km s -1 ] -1 . (9) 

The corrected X for the Ursa Major clouds is consistent with the 'canonically' adopted 
7-ray results. 
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For some clouds, however, the H 2 column density per unit CO line integral is 
significantly different from canonical values. The cloud HRK236+39 was observed in 
the 21-cm line with the Arecibo 305-m telescope, which provides the highest-resolution 
single-dish observations of an isolated, diffuse interstellar cloud ( [Reach et al. 1994]) . For 
this cloud, the linear regression of the H I, CO, and infrared maps resulted in the 'naively' 
determined value X naive = (0.16 ± 0.03). This low value is not due to a temperature effect: 
we find in the present work that the color temperature of the far-infrared excess emission 
is 16.5 K, only slightly colder than that of atomic gas. In order to explain the low X na i ve 
with a temperature difference alone, the temperature of the molecular gas would have to 
be 13.2 K. While this might seem similar to the observed 16.5 K, it is in fact clearly ruled 
out: for the same 100 /im brightness, the 240 /im emission of the cloud would be more 
than 4 times brighter than observed. For HRK236+39, comparison of the H I, CO, and 
infrared maps clearly shows that there is a component of the infrared emission that is not 
traced by H I 21-cm or CO(l — > 0) line emission. The infrared excess exists in the central 
~ 45' of the cloud, while the CO is concentrated even more to the center, in clumps with 
size ~ 2'. Similar results have been found from a comparison of infrared and CO maps of 
the Chameleon region by Boulanger et al. (1997). They found that the large clouds, with 
extensive CO emission, X naive =(0.6-1.0), which could easily be explained by a somewhat 
lower dust temperature in the molecular clouds. The large scatter in the comparison of CO 
and extinction maps — which are independent of the dust temperature — was found to be due 
to spatially coherent regions where the gas is primarily H2 but the CO abundance is low. 

EDITOR: PLACE TABLE | HERE. 

Any map of molecular gas, even over a very restricted region, would allow us to calibrate 
the infrared excess. Where it has been observed, the extinction of background starlight 
provides a temperature-independent tracer of the total column density. Comparison of 
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extinction and CO maps for high-latitude clouds yields X =(0.4-6.4) ([Hcithausen fc 
Thaddeus 1990| ; Magnani, Blitz, & Wouterloot 1988), but it is unclear whether the large 
scatter is due to real cloud-to-cloud variations or calibration difficulties in the extinction 
maps (which are obviously uncertain when the extinction is appreciably more than unity). 
The H 2 molecule is symmetric, so that its rotational transitions must be quadrupolar, 
and the lowest-energy transition has an energy of 500 K, which cannot be excited under 
typical interstellar conditions. Absorption line studies are able to trace cold, diffuse gas, 
but the paucity of suitable background stars makes comparison to large-beam observations 
problematic. A promising method is to observe the low-energy radio transitions of the CH 
molecule, because CH forms readily — but only where the H is molecular. Magnani & Onello 
(1995) have observed a set of both translucent and dark molecular clouds in the 9-cm lines 
of CH, and used the empirical correlation of CH line brightness with extinction to obtain 
the H 2 column density. They find that X =(0.3-6.8) for translucent high-latitude clouds: 
a substantial range, as was found from the extinction studies. The value of X was also 
found to vary within clouds, where it was observed for more than one position. A detailed 
comparison of our results with those of Magnani & Onello (1995) is complicated by the fact 
that the clouds are structured on scales smaller than our beam, while the CH lines have 
not been mapped over the clouds. While there is clearly an infrared excess counterpart 
of each of the clouds in their sample, only 11 lines of sight are suitable for comparison; 
we find that the ratio of 240 /zm infrared excess to H 2 column density (inferred from the 
CH line brightness) ranges from 0.4 to 2 MJy sr _1 /(10 20 cm -2 ). It is unclear whether the 
range of values is due to beam size effects or a true variation from cloud to cloud, but this 
comparison at least affords a calibration of the infrared excess independent of CO. For 
comparison, the calibration of 240 /zm emission per H 2 for the North Celestial Pole region 
yields 

4(240 fxm)/jV(H 2 ) = b/X = 1.0 ± 0.2 MJy sr" 1 (10 20 cm" 2 )" 1 (10) 
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using Table [5], which is within the range given by the CH calibration. 



EDITOR: PLACE FIGURE 12 HERE. 



4. Conclusions 

The infrared excess clouds, which were found after removing emission associated 
with the atomic interstellar medium, are predominantly molecular clouds. The dust in 
these clouds is colder than dust in diffuse atomic clouds, for the same reason that the gas 
is molecular: the outer layers of the cloud shield the center from the dust-heating and 
molecule-dissociating effects of the interstellar radiation field. These molecular clouds fill 
only a small fraction of the high-latitude sky. The column density of the infrared excess 
clouds can be estimated using the calibrations discussed above, which may be rewritten as 

JJj5,/JV(H 2 ) = (0.26 ± 0.05) x lO-^cn^MJysr -1 . (11) 

We can estimate the mass surface density of the infrared excess clouds, using the difference 
between positive and negative excess for \b\ > 20° (Fig. ^). The surface density is about 
0.3 M Q pc -2 (including 40% He by mass), roughly equally partitioned among regions of low 
and high brightness. (We adjusted for the part of the sky at |6| < 20° by assuming it has 
the same average infrared excess as the rest of the sky.) The mass density that we find is an 
order of magnitude larger than that found in an unbiased, wide-field CO survey ( Hartmami 



et al. 19981) of the northern sky, but it is comparable to a previous estimate — also based on 



CO — of the total mass of the known molecular clouds at high galactic latitude flMagnani 



et al. 1996|). We suspect that the unbiased survey missed much of the molecular gas 



because either there was no CO mixed with the diffuse H 2 or the CO rotational levels 
were subthermally excited (because of low H 2 volume density). We do not find a great 
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disparity between the northern and southern galactic hemispheres (cf. |Hartmann et al. 
|1998| ); according to the infrared excess, the southern galactic hemisphere has about 30% 
more molecular gas than the northern hemisphere. The high-latitude infrared-excess clouds 
comprise a significant fraction of the total mass of molecular gas in the solar neighborhood. 

The infrared excess maps can be used to provide a guide to regions of anomalously 
high extinction. The extinction calculated from the H I column density (cf. [Burstein and| 
Heiles 19"78| ) can significantly underestimate the total extinction in these regions. If we 
assume that the dust in the infrared-excess clouds has the same extinction cross-section per 
H-nucleus as diffuse clouds for which this quantity has been measured (Bohlin, Savage, & 
Drake 1978), then the extinction due to infrared excess clouds can be obtained from Fig. 0, 
after multiplying the 100 /xm surface brightness by Ay/I^ Q ~ 0.027 mag MJy -1 sr. 

WTR was supported by NASA's COBE project during the early part of this work, 
through a contract to the Universities Space Research Association. In searching for 
counterparts for the objects found in this work, we use the Sky View facility at NASA 
Goddard Space Flight Center; Sky View is supported by NASA ADP grant NAS 5-32068. 
We also made use of the SIMBAD database, which is maintained by the Centre de Donnees 
astronomiques de Strasbourg (CDS), France. We thank Xavier Desert, Frangois Boulanger, 
Rick Arendt, Janet Weiland, Tom Sodroski, and Mike Hauser for their comments and 
suggestions. 
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Fig. 1. — Infrared surface brightness at 240 /zm observed by DIRBE plotted versus the atomic 
gas column density from the Leiden-Dwingeloo 21-cm line survey for the Ursa Major/Ursa 
Minor/ Camelopardalis region. In the top panel, nearly the full dynamic range of H I and 
infrared brightness is shown; 56 points, out of 2393 total, with N(R I) > 1.6 x 10 21 cm -2 
were excluded. Error bars are shown for a few points for illustration. The vertical dashed 
line indicates the threshold column density below which a linear fit was performed. The 
bottom panel shows the same data, but with restricted dynamic range covering only the 
points used in the fit, iV(H I) < 3.75 x 10 20 cm~ 2 . The linear fit to the low-column density 
points is shown in both panels as a diagonal solid line. In this paper we concentrate on the 
positive residuals from the linear fits, which we attribute to dust associated with ionized and 
molecular gas. 

Fig. 2. — All-sky maps of the slope (a) and offset (b) of the correlation between 100 jum 
surface brightness and H I column density. The maps are shown in a galactic-coordinate 
Mollweide projection with the galactic center at the center. Meridians are drawn every 60°, 
and parallels are drawn every 30°. The greyscale ranges are labeled to the right of the images; 
for panel (a) in MJy sr _1 / (10 20 cm -2 ), and for panel (b) in MJy sr _1 . 

Fig. 3. — Histograms of the slopes of the infrared-H I correlations for sky patches with 10° 
radius. Each panel shows the results from a different DIRBE waveband: (a) 100 /im, (b) 140 
/zm, (c) 240 /im. The 'outlying' values of the slope, defined as being more than 1.5 times 
the rms away from the mean, are shaded. 

Fig. 4. — Similar to previous figure, but this time for the offsets of the infrared-H I 
correlations for sky patches with 10° radius. Each panel shows the results from a different 
DIRBE waveband: (a) 100 /im, (b) 140 fj,m, (c) 240 /im. The sky patches with 'outlying' 
slope are shaded. 
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Fig. 5. — All-sky map of the residual 240 fim surface brightness after subtracting the emission 
correlated with the H I column density. The map is shown in a galactic-coordinate Mollweide 
projection as Figure 2; meridians are drawn every 60° and parallels every 30°. The greyscale 
range is illustrated by the bar, with units MJy sr^ 1 , to the right of the map. 

Fig. 6. — Histograms of the infrared excess at 100 /xm. The upper histogram shows the 
distribution of the positive values of infrared excess, and the curve beneath it shows the 
distribution of the negative values. The positive values clearly dominate above an infrared 
excess of 1 MJy sr _1 . 

Fig. 7. — An atlas of the 100 fim infrared excess in galactic-coordinate Cartesian projections, 
for absolute galactic latitudes 20° to 50°. Each map covers 120° of longitude, with the centers 
at longitudes 0° (a), 120° (b), and 240° (c). The upper part of each page shows positive 
galactic latitudes, and the lower part shows negative galactic latitudes. In each panel, the 
most significant clouds are identified and labeled using the catalog of Magnani et al. (1997) 
or, for new clouds the 'DIRlllibb' cloud name. The greyscale ranges from -0.2 (white) to 
5 (black) MJy sr _1 at 100 /im wavelength. Contours are drawn at 0.3, 1, 2, 3, and 4 MJy 
sr _1 (but the highest levels are not present in every map). The large black areas are masked, 
either because of lack of H I data, presence of the Magellanic clouds, or presence of extremely 
bright and structured infrared emission. 

Fig. 8. — Maps of the 100 fxm infrared excess at the galactic poles, in orthographic 
projections. Panel (a) covers the north galactic pole and panel (b) covers the south galactic 
pole. As in the previous figure, the most significant infrared excess clouds are labeled. The 
greyscale range, from white to black, is -0.2 to 3 MJy sr _1 . Contours are drawn at 0.1, 1, 
0.5, 1, and 2 MJy sr _1 . (The contours are at lower values in this figure than in Figure 7, in 
order to bring out fainter features near the poles.) The large black stripe in panel (b) is due 
to a region where the H I 21-cm data are poorly sampled. 
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Fig. 9. — Surface brightness map at 60 /im of the warm infrared excess cloud DIR015+54, 
made from the IRAS Sky Survey Atlas. The minimum surface brightness in the map is -1.83 
MJy sr" 1 (negative because of oversubtraction of zodiacal light), and the contours are at 
-1.35, -1, and -0.65 MJy sr" 1 . 

Fig. 10. — Histograms of color temperatures of MBM clouds (top) and DIRBE clouds 
(bottom) obtained from the infrared excess at 100 and 240 fxm. The color temperature 
of emission associated with the atomic gas is shown as a dashed line. Both the MBM clouds, 
which are known to be molecular, and the new infrared excess clouds are colder than atomic 
clouds, most likely because they are optically thick to the high-energy part of the interstellar 
radiation field. 

Fig. 11. — Comparison of the infrared excess (smooth contours) and CO(l — > 0) line integral 
for the Ursa Major-Ursa Minor-Camelopardalis region. Regions where the CO was observed 
are surrounded by thick lines and filled with thin diagonal lines. The two maps are in very 
good agreement where data for both maps exist. The infrared excess predicts, accurately, 
the location and extent of the molecular clouds. Furthermore, there are no 'extra' infrared 
excess clouds predicted in the regions where CO was observed but not detected. 

Fig. 12. — Correlation of the infrared excess, /fgo (MJy sr -1 ), and CO(l — > 0) line integral, 
W(CO) (K km s _1 for the North Celestial Pole region. Each point corresponds to an 
independent 40' x 40' pixel. 
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Table 1. Infrared to H I Correlations for Two Regions (or Cells) 



(/,6) = (140°,35°) (;,6)=(40°,35°) 
Wavelength Slope Offset \l Slope Offset xl 

[MJy sr- 1 / 10 20 cm- 2 ] [MJy sr" 1 ] [MJy sr" 1 / 10 20 cm- 2 ] [MJy sr" 1 ] 



100 0.56 ±0.01 0.63 ± 0.02 1.25 0.89 ± 0.01 0.33 ± 0.04 1.77 

140 ^m 1.12 ±0.06 1.01 ±0.14 1.06 1.45 ±0.06 0.44 ±0.21 0.91 

240 /^m 1.02 ±0.03 0.62 ± 0.08 1.26 1.07 ±0.04 0.42 ±0.12 0.91 



Table 2. Statistics of Slopes and Intercepts 3, 



Wavelength 


Slope 


Offset 


Offset w/o outliers b 




[MJy sr" 1 / 10 20 cm" 2 ] 


[MJy sr" 1 ] 


[MJy sr" 1 ] 


100 ^ni 


0.62 ±0.08 


0.76 ±0.37 


0.71 ±0.32 


140 fim 


0.95 ±0.17 


1.11 ±0.63 


1.10 ±0.52 


240 ^ni 


0.72 ±0.11 


0.87 ±0.55 


0.83 ±0.41 



a mean and rms dispersion, from cell to cell 

b excluding those regions for which the slope deviates by more than 1.5 times 
its dispersion from its mean value 
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Table 3. Infrared Excess from Known High-Latitude Molecular Clouds 



IN d-lllt, 






S1ZG 




zpex b 
r 100 










[Kl 


[k Tvl 


DC001.4-21.6 


1.4 


-21.6 


0.6 


16.2 


1.31 


R Cor Bor 


4.0 


-25.5 


2.0 


17.1 




LDN 134 


5.0 


36.0 


3.0 


14.8 




MBM 45 


9.8 


-28.0 


1.0 


14.2 


0.72 


MBM 39 


11.4 


36.2 


1.9 


16.6 


1.56 


MBM 40 


37.6 


44.7 


2.2 


15.9 


0.48 


MBM 46/48 


40.3 


-35.3 


2.0 


16.9 


2.12 


HD 210121 


57.5 


-44.0 


4.0 


16.3 


3.57 


G61-34 


61.3 


-34.0 


1.3 


15.2 


1.16 


MBM 49 


64.7 


-26.8 


1.0 






MBM 50 


70.0 


-30.7 


5.5 


21.0 


1.85 


G72+25 


71.0 


25.5 


1.0 


17.7 




MBM 51/52 


74.0 


-51.3 


4.0 


17.0 


0.66 


MBM 55 


89.2 


-40.9 


3.4 


14.6 


1.77 


MBM 41 /44 


90.0 


38.9 


1.8 


15.7 




MBM 53 


93.0 


-32.2 


4.4 


19.0 




MBM 54 


93.0 


-37.5 


4.3 


15.9 




HRK 94.8+38 


94.8 


37.5 


1.2 






MBM 56 


102.7 


-28.1 


2.5 


17.2 


1.24 


IVlrSlVl 1 


lino 
1W.Z 


-41. Z 


U.O 


17 Q 
Li .O 


U.O / 


MBM 2 


117.7 


-52.6 


1.9 


19.0 


1.28 


HSVMT 1 


119.1 


28.5 


1.8 


15.9 


0.61 


Polaris 


125.0 


28.0 


5.0 


15.1 


0.79 


HSVMT 12/14 


125.5 


32.5 


1.0 


16.0 


0.33 


MBM 3 


130.0 


-46.8 


2.5 


15.4 


1.21 


MBM 4 


134.2 


-45.0 


2.0 


14.7 


0.63 


HRK 135+51 


135.6 


51.3 


1.5 


16.2 


0.62 


HRK 135-69 


136.4 


-68.3 


1.8 


17.9 


2.32 



Table 3 — Continued 



Name a 


7 

t 


7 




size 
[degj 


1 (240/11)0) 
[K] 


^100 

n t l 

[kJyJ 


HRK 135+54 


136.5 


54.8 


1.1 


17.5 


2.86 


HRK 140+48 


140.5 


48.0 


0.7 


24.4 


0.74 


MBM 27/29 


142.0 


35.0 


2.2 


15.6 




MBM 30 


142.2 


38.2 


3.4 


15.2 


0.50 


Camelopaxdalis 


143.5 


24.0 


2.9 


16.0 


0.61 


MBM 5 


145.0 


-49.9 


3.8 


13.8 


0.60 


MBM 6 


147.0 


-39.0 


1.9 


15.0 


1.07 


MBM 31/32 


147.0 


40.0 


4.2 


15.9 




MBM 7-8 


150.9 


-38.3 


1.5 


14.8 


2.87 


HSVMT 27 


153.6 


36.9 


2.8 


15.7 


2.20 


G 155-40 


154.7 


-39.8 


1.4 


15.2 


0.91 


MBM 26 


156.4 


32.6 


2.0 


18.1 




MBM 11/14 


159.0 


-33.6 


3.5 


15.3 


0.62 


MBM 17 


165.7 


-25.6 


0.5 






MBM 16 


171.2 


-37.7 


0.5 


15.5 


0.77 


MBM 23/24 


172.0 


26.9 


0.8 


15.5 


0.84 


MBM 25 


173.6 


31.2 


2.0 


15.9 




MBM 19 


186.0 


-29.9 


0.5 






MBM 18 


189.1 


-36.0 


0.5 






a ill ti it -1 r 

MBM 15 


191.7 


-51.3 


3.5 


17.7 




HRK 192-67 


192.0 


-67.5 


1.4 


18.5 


2.29 


CB 28 


204.0 


-25.2 


1.5 


14.9 


0.28 


MBM 21-22 


208.3 


-28.0 


2.0 


19.6 


0.53 


G211+63 


210.8 


63.1 


2.0 


18.4 




MBM 20 


210.9 


-36.5 


2.5 


15.4 


2.27 


HRK 225-66 


227.0 


-66.1 


2.5 


15.6 


2.69 


HRK 228-29 


229.0 


-28.5 


1.7 


16.4 




HRK 236+39 


236.0 


38.0 


4.0 


16.5 


0.22 
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Table 3 — Continued 



Name 3, 


I 


b 


size 


T(240/100) b 


TP ex b 
^100 








[deg] 


[K] 


[kJy] 


KM 273+29 


272.5 


29.7 


2.0 


17.9 


0.65 


KM 293-31 


291.0 


-31.0 


2.2 






KM 300-24 


301.1 


-24.5 


2.0 






MBM 33 


359.1 


36.7 


1.2 


15.7 





a Cloud names mostly from the compilation of high-latitude 



molecular clouds (Magnani et al. 1996). MBM: extinction 



patches on Palomar Observatory Sky Survey (Magnani et al 



1985), KM, DC: extinction patches in European Southern 



Observatory sky survey (Keto and Myers 1986; Hartley et al 



1986), HRK: infrared cirrus patches in IRAS 100 /im maps 



flHc-iles et al. 1988[ ) 

b Values listed as in clouds for which the COBE/DIRBE 
data indicates strong, unresolved structure (consistent with a 
bright point source.) 
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Table 4. Catalog of Unidentified Infrared Excess Clouds 



Name I b size T(240/100) Ffg 

[deg] [K] [kJy] 



DIR002+31 


1.6 


-30.8 


1.2 


18.4 


0.89 


DIR009+30 


9.1 


-30.1 


2.3 


16.5 


2.50 


DIR013+40 


12.8 


40.0 


1.5 


17.1 


1.31 


DIR018+37 


17.8 


-36.5 


2.5 


12.2 


0.04 


DIR020+45 


19.5 


-44.5 


1.7 


18.9 


0.69 


DIR021+52 


21.0 


52.0 


1.5 


17.5 


0.29 


DIR025+35 


25.0 


35.0 


3.5 


21.8 


2.96 


DIR027+31 


26.5 


-30.5 


1.0 


15.3 


0.52 


DIR028+54 


28.3 


53.7 


1.0 


17.3 


0.29 


DIR029+25 


28.5 


25.0 


2.6 


16.6 


2.15 


DIR029+30 


28.5 


30.0 


1.1 


16.7 


0.72 


DIR034+26 


34.0 


26.0 


3.0 


18.6 


3.21 


DIR046+37 


45.5 


-36.5 


1.6 


16.1 


1.24 


DIR046+33 


45.7 


-33.0 


2.2 


22.4 


0.72 


DIR048+25 


47.7 


24.5 


2.1 


16.9 


1.49 


DIR048+38 


48.0 


38.0 


2.5 


18.0 


1.56 


DIR060+26 


59.7 


-26.0 


1.7 


22.3 


1.07 


DIR061+22 


61.0 


22.0 


4.0 


26.1 


3.57 


DIR070+23 


69.7 


22.5 


1.2 


17.2 


0.63 


DIR071+43 


71.0 


-42.5 


3.0 


15.5 


3.14 


DIR072+34 


71.5 


-34.0 


1.4 


15.9 


0.43 


DIR077+37 


76.5 


-37.0 


3.3 


18.5 


1.75 


DIR081+39 


81.3 


38.5 


2.5 


16.2 


0.80 


DIR087+29 


86.5 


28.5 


3.5 


23.1 


1.07 


DIR096+23 


96.2 


23.0 


2.0 


19.3 


1.17 


DIR098+44 


98.0 


-44.0 


2.5 


17.4 


0.97 


DIR105-31 


105.0 


-31.0 


1.8 


15.3 


1.16 


DIR105-38 


105.0 


-38.0 


2.7 


15.7 


2.48 
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Table 4 — Continued 



Name I b size T(240/100) Ffg 

[deg] [K] [kJy] 



DIR108-27 


108. 


.0 


27. 


.0 


2 


.7 


18. 


1 


2.12 


DIR1 17-44 


116. 


.5 


-44 


.0 


2 


.0 


18. 


.0 


0.70 


DIR120-28 


120. 


.0 


-28 


,0 


2 


.4 


17. 


.7 


1.28 


DIR121-45 


120. 


.9 


-45 


.0 


1 


.8 


15. 


.5 


0.57 


DIR126-37 


126. 


.0 


37. 


.0 





.9 


16. 


,7 


0.30 


DIR132-30 


131 


5 


-29 


.5 


1 


.9 


18. 


.4 


0.61 


DIR134-36 


134. 


.0 


-35. 


.5 


1 


.5 


16. 


.3 


0.62 


DIR135-38 


134. 


.5 


38. 


.0 


1 


.4 


23. 


.4 


0.28 


DIR135-41 


134. 


.7 


40. 


.6 


1 


.4 


17, 


.0 


0.34 


DIR140-45 


140. 


.0 


-44. 


.5 


2 


.0 


19. 


.8 


0.63 


DIR150-29 


150. 


.2 


-28. 


.5 


1 


.7 


16. 


2 


0.89 


DIR152-47 


152. 


.0 


-46. 


.5 


3 


.0 


16. 


2 


2.74 


DIR164-44 


163. 


.5 


-44. 


.0 


1 


.6 


17. 





0.61 


DIR1 72-42 


172. 


.0 


-41. 


.5 


2 


.3 


15. 


6 


2.87 


DIR177-33 


177. 


.2 


33 


.0 


1 


.5 


16. 


.3 


0.56 


DIR1 79-49 


179. 


.0 


-49 


.0 


3 


.5 


20. 


.4 


1.68 


D1R184-26 


183 


.5 


26. 


.0 


1 


.8 


17, 


,8 


0.66 


DIR187-43 


187. 


.0 


-42 


,5 


1 


.7 


20. 


.4 


0.58 


DIR196-24 


196. 


.0 


24. 


,2 


1 


.0 


17. 


.2 


0.46 


DIR198-32 


198. 


.0 


32 


.0 


1 


.5 


15. 


,8 


0.61 


DIR201-24 


200. 


.5 


23 


.6 


1 


.3 


19, 


9 


0.39 


DIR203-32 


202. 


.5 


-31 


.5 


2 


.0 


14. 


9 


0.33 


DIR204-37 


203. 


.5 


-36 


.5 


1 


.8 


15. 


.6 


0.79 


DIR216+27 


216. 


,3 


26. 


.5 


2 


.2 


16. 


2 


0.34 


DIR223+37 


223. 


.0 


-37 


.0 


2 


.7 


22. 


3 


1.21 


DIR234+37 


234. 


.0 


-37 


.0 


1 


.1 


21. 


2 


0.65 


DIR237+44 


237. 


.0 


-44 


.0 


1 


.2 


19. 


,7 


0.72 


DIR239-25 


239. 


.0 


-25. 


.0 


4 


.0 


20. 


5 


3.34 



Table 4 — Continued 



Name I b size T(240/100) Ffifo 

[deg] [K] [kJy] 
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ID. 
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.Z 


1 Qfi 

i.yo 


PUT? 97/1 i/IR 
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Z ( 4. 


.u 


-40 


Q 
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Table 5. Multiple Linear Regression between Infrared, H I, and CO Maps 



Wavelength a = I V /N(H I) b = I v /W(CO) xl X = N(H 2 )/W(CO) 

[MJy sr- 1 / !0 20 cm- 2 ] [MJy sr" 1 / (K km s" 1 )] [10 20 cm- 2 /(K km s" 1 )] 

naive true 



100 Mm 0.48 ±0.02 0.34 ±0.03 1.1 0.35 ± 0.05 1.3 ±0.3 

140 /im 1.06 ±0.08 1.24 ±0.12 0.8 0.58 ±0.11 1.5 ±0.4 

240 Mm 0.94 ±0.06 1.33 ±0.10 1.2 0.71 ± 0.11 1.3 ±0.2 



a For the Ursa Major/Ursa Minor/Camelopardalis region 
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